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A  THERMODYNAMIC  MODEL  FOR  C02  LASERS 
by  D.  C.  Dumitras 

Central  Physics  Institute,  Bucharest  P.O.  Box  5206 
Inst,  for  Physics  and  Technology  of  Radiation  Equipment 

1.  Introduction 

In  ord  r  to  optimize  the  functioning  of  a  CC>2  laser 
with  wave  guide,  it  is  important  to  be  able  to  calculate 
the  dependence  of  the  small  signal  gain  as  a  function  of 
a  number  of  parameters  such  as:  gas  mixture,  total  pressure, 
discharge  current,  tube  diameter,  temperature  of  the  cooling 
jacket,  etc.  If  the  kinetic  temperature  and  the  population 
of  the  laser  levels  are  known  as  a  function  of  these  para¬ 
meters,  the  gain  equations  are  uniquely  determined.  The 
complexity  of  the  energy  level  diagram  for  the  C02  molecule, 
the  large  number  of  rates  of  colllsional  excitation  and 
relaxation  and  their  dependence  on  temperature,  diffusion 
and  relaxation  processes  through  collisions  with  the  walls 
of  the  enclosure,  the  great  number  of  additives  that  parti¬ 
cipate  in  the  population  inversion,  add  to  the  difficulty 
of  the  analytical  treatment  of  the  laser  effect  and  of  the 
gain,  saturation  and  oscillation  bandwidth. 

Since  its  inception,  a  number  of  models  were  proposed 
based  on  the  rate  equations  relating  the  vibrational  modes 
of  the  C02  and  N2  molecules.  Thus,  Moore  et  al  [1]  have  di¬ 
vided  the  energy  levels  in  three  groups  (each  characterized 
by  a  single  vibrational  temperature),  and  solved  the  system 
of  equations  for  the  continuous  wave  regime.  A  comparison 
with  experiment  shows  that  this  is  a  good  simplification. 
Previous  works  considered  separately  the  vibrational  levels 
of  nitrogen  [2],  [3],  and  4  or  5  groups  including  also  the 
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vibrational  levels  of  Co[4], 

The  new  elements  in  this  work  are:  small  tube  diameters, 
high  pressures  and  current  densities  characteristic  to  the 
wave  guide  lasers.  These  differences  imply  faster  excitation 
and  relaxation  rates  and  also  collisional  widening  of  emis¬ 
sion  lines.  The  most  recent  determinations  of  the  excitation 
and  relaxation  rates  and  their  temperature  dependence  are 
utilized.  This  data  was  not  available  when  [1]  -  [4]  were 
published. 

The  main  limitation  of  works  dealing  with  thermodynamic 
models  for  C02  lasers  are  using  calculated  excitation  and  relax¬ 
ation  rates  [2-5],  which  can  differ  appreciably  from  the  measu¬ 
red  ones;  the  gas  temperature  is  determined  from  the  wrong  relat¬ 
ionship  [2-5];[l]  neglects  relaxation  through  diffusion  and  col¬ 
lision  at  the  walls  of  the  enclosure  and  is  limited  to  a  3 
level  system;  none  of  these  works  take  into  account  the  gas 
pressure  dependence  of  the  diffusional  relaxation  rate,  and 
neglect  the  losses  inherent;  when  guiding  radiation  through 
small  diameter  tubes;  [5]  is  limited  to  a  fixed  excitation  rate 
and  only  one  diameter  of  the  discharge  tube  and  does  not  take 
into  account  the  saturation  processes.  Some  recent  attempts 
[63],  [7]  eliminate  some  of  these  problems  but  the  system  of 
equations  used  is  still  simplistic  (3  level  system). 

The  model  presented  here  is  used  for  lasers  with  wave  guide. 
It  can  be  used  also  for  other  C02  lasers  with  axial  discharge 
and  gas  flow.  For  lasers  without  gas  flow,  the  modifications 
of  the  rate  equations  are  due  to  the  kinetic  temperature  of  the 
gas  and  eventual  gaseous  components  which  participate  in  the 
inversion  in  the  C02  molecule  (especially  the  carbon  monoxide). 

2.  Physical  Model 

The  lower  vibrational  energy  levels  for  the  C02  and  N2 
system  are  shown  in  figure  1. 
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These  levels  are  divided  into  four  groups:  group  0, 

including  the  fundamental  levels  (unexcited) of  the  C0?  and  N5 

0  d  ^ 

molecules,  with  populations  and  respectively;  group  1, 
including  the  C02  molecules  excited  in  the  modes  of  symmetric 
elongation  and  deformation,  with  population  N^;  group  2,  C02 
molecules  with  asymmetric  elongation  modes  excited,  population 
N2;  group  3,  containing  nitrogen  molecules  with  one  or  more 
of  the  vibrational  levels  excited  and  population  N^.  The  total 
densities  of  C02  and  nitrogen  molecules  are  and  respectively 


Dividing  the  system  C02  and  N2  in  four  groups  seems  just¬ 
ified  by  the  rapid  relaxation  of  the  vibrational  energy  within 
each  group,  at  a  rate  much  larger  than  that  of  exchange  of  energy 

between  groups.  Thus  C0~  molecules  in  state  00°2  relax  to 

d.  6  l  1 

state  00°1  at  a  rate  larger  than  5  x  10  s-  torr-  [8],  and 
the  nitrogen  molecules  relax  from  the  upper  levels  to  state 
v  =  l  in  approximately  lps  [9].  Also  the  relaxation  of  the  vi¬ 
brational  energy  between  levels  in  group  1  is  rapid,  it  is  thus 
possible  to  characterize  the  whole  group  by  a  single  vibrational 
temperature.  Groups  1-3  are  excited  through  electronic  colli¬ 
sion  with  rates  R-^  -  R^  and  relax  through  diffusion  and  colli¬ 
sions  with  the  walls  of  the  enclosure  with  rates  k.  -  k0  . 

IP  3p 

Group  1  is  excited  and  relaxes  collisionally  with  rates  k^Q 
and  kQ1  respectively,  and  group  3  relaxes  through  electronic 
collisions  at  a  rate  k^g.  Vibrational  energy  is  exchanged  be¬ 
tween  groups  3  and  2  through  a  rapid  V-V  transfer,  character¬ 
ized  by  a  rate  k^  and  the  inverse  rate  k^,  and  between  groups 
2  and  1  through  collisions  with  molecules  and  atoms  (rates  k21 
and  k12),  collisions  with  discharge  electrons  (rate  k2„),  through 
induced  emission  (rate  W21)  and  induced  absorption  (rate  W^2). 

The  inverse  rate  k^  between  groups  2  and  3  is  taken  to  be  equal 
to  the  direct  rate  k^2*  Rate  k^0  (due  to  collisions  with 
molecules,  atoms  and  electrons)  is  neglected,  because  such  a 
collision  is  much  more  likely  to  relax  a  C02  molecule  from 
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Figure  1.  Diagram  of  lower  vibrational  energy  levels  for 
the  C0„  and  system,  grouped  for  the  4  level 
approximation. 

1  -  energy; 


state  00°l  to  one  of  the  states  in  group  1,  rather  than  the 
fundamental  level  [9].  We  assume  also  that  the  degree  of 
dissociation  of  the  CC>2  molecule  is  small  (justified  for  the 
case  of  lasers  with  gas  circulation). 

2.1.  Kinetic  Temperature  of  the  Gas 

To  calculate  the  gas  temperature  in  the  discharge  one 
needs  to  solve  the  Fourier  heat  equation  for  the  continuous 
wave  generation  regime 

«  4  f  -  </;  J-uii,  —  0, 

(1) 

where  T  is  the  gas  temperature  in  K,  g  is  the  rate  of  change 

U 
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of  the  heat  due  to  sources  and  radiators  expressed  in  Wm-^ 
and  mix  is  the  thermal  conductivity  of  the  gas  mixture  in 
WK_1rrf 1 . 


The  source  of  heat  is  that  part  of  the  electron  energy 
which  is  transferred  to  the  translational  energy  of  the  mole¬ 
cules.  Thus  the  radial  distribution  of  the  heat  source  is 
the  same  as  the  radial  distribution  of  the  electrons  in  the 
discharge . 


From  the  theory  of  the  positive  column  the  radial  distri¬ 
bution  of  the  electron  density  is  a  zero  order  Bessel  function, 
if  the  distribution  is  not  affected  by  the  strong  radial  temp¬ 
erature  gradient 


Solving  the  heat  equation  with  this  q  in  cylindrical  coordinates 
with  boundary  condition  T(a)  =  T  ,  for  a  cylinder  of  radius  a 
and  length  L  (length  of  discharge),  one  obtains 


T 


'•'HIV 


(2,405  r!«). 


(3) 


If  we  also  consider  the  temperature  gradient  in  the  wall  of 
the  discharge  tube  (exterior  radius  a')  then  [7] 


T 


,  r  i  ,  /«'\i 

>■ .  -»*.*(.)><*•*“'«• 


where  K  is  the  thermal  conductivity  of  the  wall  material. 
In  the  center  of  the  tube  (r=0)  the  gas  temperature  can  be 
calculated  from 


r 


-  t. 


0,1729 

ini  lv 


(5) 


The  total  amount  of  heat  which  increases  the  translational 
energy  of  the  molecules  per  unit  time,  Q,  is  obtained  by  inte 
grating  the  source  q  over  the  volume 


■MjjpgMHMiMI 


V  -  \  </o')uC-i,  »Ur»  r/«)  ill'  =  1,3.-. ttl*  q0I*a*. 


(6) 


Thus 


T  =  T 


P 


V  r»42M  r  o,:mwg  ,j«:y 

/.  .  irinx  ip  i  ^  . 


(7) 


(Works  [2]  -  [5]  have  used  this  formula,  'out  the  proportionality 

coefficient  for  T  -  T  was  2.4  times  smaller). 

]? 


In  calculating  the  temperature  of  the  gas  mixture  in  the 
center  of  the  tube,  one  has  to  take  into  account  that  not  all 
the  input  power  is  turned  into  heat:  part  is  turned  into  laser 
power;  part  to  excite  the  molecules  which  in  turn  through  dif¬ 
fusion  and  collisions  lose  it  to  the  walls;  and,  in  the  case 
of  lasers  with  gas  circulation,  part  of  the  excited  molecules 
leave  the  discharge  zone,  before  their  energy  is  turned  into 
heat.  Because  of  these,  equations  (5),  (7)  have  to  be  adjusted 
by  introducing  a  conversion  coefficient 


r-  t, 

T-  T, 


<•.  1 7bi> 

'•IIIIV 


v  nudTt 


(3) 

(9) 


For  closed  C02  lasers,  >)**  varies  between  0.5  and  0.B  [10], 
and  for  those  with  gas  circulation  between  0.05  and  0.3. 


One  can  see  from  (9)  that  T  -  T  is  proportional  with  the 

* 

input  electrical  power  per  unit  length  (Q/L),  thus  the  axial 
temperature  increases  with  increasing  current  and  pressure 
and  for  same  Q/L  is  independent  of  the  radius  of  the  tube. 


Laderman  and  Byron  [11]  found  the  exact  solution  for  the 
heat  equations,  and  the  diffusion  equation  for  the  electrons. 
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They  showed  that  the  radial  distribution  of  the  electron 
density  differs  only  slightly  from  the  Bessel  profile  used  here 
for  temperatures  at  the  center  of  the  tube  less  than  600  K.  They 
also  found  the  temperature  at  the  center  of  the  tube  to  be 
a  function  of  the  input  power  per  unit  length  only. 

The  sources  of  heat  are  going  to  be  estimated  in  two 

ways . 


First  by  using  eq.  (3)  assuming  that  the  gas  temperature 
is  due  to  the  energy  transferred  from  the  electrons  to  the 
molecular  vibrations  and  subsequent  relaxation  to  ground  level. 
From  fig.  1  one  obtains 

</„  —  l(/»\  -  ) A n  E3S"  —  t.V3J  A3,  (10) 


where  Ne  is  the  electron  density,  E 

—  20 

of  quantum  ,j3  ,  E^  =  4.6653  x  10 


3 

T 


is  the  vibrational  energy 


Second,  using  eq.  (9)  and  the  input  power  per  unit  length 
estimated  from 


L 


IE  —  IS- 


E 

S 


E 

T 


(ii) 


where  N  and  p  are  the  total  density,  and  pressure  of  the  gas 
mixture,  E  is  the  axial  electric  field  and  I  the  discharge 
current. 


Lowke  eb  al  [12]  have  solved  numerically  the  Boltzmann 
equation  and  derived  the  transport  coefficients  for  the  elec¬ 
trons  in  the  gas  mixtures  specific  to  C0?  lasers.  For  j/p  ® 

_  "5  —2  —1  ^ 

10  J  -  10  A  cm  torr  typical  for  C0o  lasers  with  wave  guide 

—  T  6 

(j  is  the  current  density)  E/N  varies  between  (2-5)  x  10“  V 

n  _  1  &  o 

cm  .  In  the  following  we  assume  E/N  =  3  x  10“  V  cm  thus  (11) 
becomes 

V  I-  ('V  tn)  i  t  ( m A )./»  ( torr). 


(12) 
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is  approx 


The  thermal  conductivity  of  the  gas  mixture,  , 

imated  using  a  formula  given  'ey  has  on  and  Saxena  [13] 


K  \ 

hi r# 


(13) 


where  n  is  the  number  of  gases  in  the  mixture,  and  >, 

are  the  thermal  conductivity  and  molar  fraction  of  the  i-th 
component,  G . ,  is  a  function  of  the  molecular  weight  of  comp- 

iiC 

onent  i  and  the  viscosities  of  components  i  and  k.  For  the 
CQp  +  N?  +  He,  and  the  mixture  [3] 

A.  , ...  _ Inj _ 

1  D.-'* I  VSj/  l  1G  'Ji !»._  Vn.  "  >n. .  y  v  ^  ,  ■  . 

_ All.-  _  ■  ^  ' 

1  -  •!,  1  yin .  'jut  -•<  Vnj/Vh>- 

where  ,  is,  ,  >.n,  are  the  conductivities  of  the  pure 

gases  and  >iut  ,  ,  yn,  are  the  molar  fractions. 


The  thermal  conductivities  of  the  pure  gases  (in  WK  xm  “) 
heir  temperature  dependence  were  measured  cy  reeth  and 
:  >,n:  -  <..x-  to-*  r >.N,  ~ •  nr'r,ai  xh..  -- 


:>.m  •  k,-:i  Ti7' 


The  thermal  conductivity  of  the  wall  ^p  ,  is  for  Pyrex 
glass  (borosidicate  732  )  '>•„  =  0.172  T  0,33  [15  ]?  for  alumi¬ 

na  -  37.7  [16]  and  for  beriliia  =  219.3. 


The  temperature  dependence  of  the  thermal  conductivity 
of  the  mixture  CO2  +  ^  and  He  is  shown  in  fig.  2.  Figure  3. 
presents  the  thermal  conductivity  of  CC2  +  He  and  + 

He  mixtures  as  a  function  of  the  molar  fraction  of  the  He.  Du 
to  the  high  thermal  conductivity  of  He...  the  conductivity  of 
the  mixture  is  higher  the  higher  the  fraction  of  He  and  the 
higher  the  temperature. 


Due  to  the  temperature  dependence  of  the  conductivity 
of  the  mixture,  the  temperature  is  calculated  from  equation 
(8)  and  (9)  by  iteration.  We  assume  initially  T  =  Td ,  and 
then  >.„,v  and  the  temperature  of  the  mixture  are  calculated. 
Then  a,„u  is  recalculated  for  the  new  temperature  and  the 
process  is  repeated.  The  iteration  is  continued  until  the 
difference  between  two  successive  calculated  temperatures  is 
less  than  1  K. 


Figure  4  shows  the  difference  between  the  temperature 
at  the  center  of  the  tube  and  the  wall  temperature  as  a  function 
tf  the  input  power  per  unit  length  for  the  tXh  :  He  -  l :  ,j ;  kj  CO,  : 

:  IU>  -  i  :  i,i)«:  4  mixtures.  The  continuous  lines  are  cocaine o 
-sing  (if),  the  cashed  lines  are  due  to  Ladermai  and  Byron  (11] 
for  the  same  mixtures.  The  circles  are  experimental  results 
by  the  same  authors  for  the  i  :  it*  mixture.  The  agreement 

between  our  calculation  and  the  experimental  data  is  excellent. 


T  -  Tp  is  a  strong  function  of  the  wall  material,  as  can 
be  seen  from  figure  5.  Thus,  the  curve  a' =4  mm  is  for  Pyrex 
glass  and  the  curve  a'  =  1  mm  is  for  all  purposes  valid  for 
alumina  or  berillia  tubes  (or  when  the  temperature  gradient 
in  the  wall  is  neglected).  For  a  Pyrex  glass  tube  with  a  = 

1  mm  and  a*  =  4mm  and  a  mixture  CO^ :  W2 :  He  =  1:1:3,  tbe  temp¬ 
erature  graident  in  the  wall  for  T  -  T  -  500  K  is  over  20" 

tm 

of  T  -  T  . 

P 
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Figure  6  shows  the  difference  between  the  temperature  at  the 
center  of  the  tube  and  the  wall  temperature  as  a  function  of  the 
total  pressure,  for  various  discharge  currents.  The  continuous  lines 
are  calculated  from  (9),  circles  are  calculated  from  (8),  for  the 
same  discharge  conditions.  The  populations  in  (8)  are  obtained  by 
solving  on  the  computer  the  system  of  equations  presented  in  section 
3.  The  ratio  n#/n**  was  taken  to  be  1.45  for  1=2  mA,  1.73  for  1=5  mA 
and  1.97  for  1=8  mA.  The  results  obtained  via  these  two  methods  are  in 
very  good  agreement.  In  the  following  the  temperature  of  the  gas 
mixture  at  the  center  of  the  tube  is  calculated  from  (9)  only,  due  to 
the  fact  that  the  computer  time  is  much  less. 


Figure  5:  Difference  between 
temperature  at  the  center  of 
the  tube  and  at  the  wall  as 
a  function  of  input  power  per 
unit  length  with  and  without 
temperature  gradient  in  the  wall. 


Figure  4.  Difference  between 
temperature  at  the  center  of 
the  tube  and  at  the  wall  as 
a  function  of  input  power  per 
unit  length  for  the  mixtures 
CC>2 :  N2  :  He  =  1:3:16  and  1:1: 

66:4. 


Figures  7-9  show  the  temperature  of  the  gas  mixture  at  the 
center  of  the  discharge  tube  for  0  **  --  <M.v ,  as  a  function  of  the 
total  pressure,  discharge  current  and  respectively  the  ^^!all ' 
temperature  (temperature  of  the  cooling  agent).  Figures  10-11 
show  the  calculated  temperature  versus  the  molar  fraction  of  N2 
and  He  respectively.  The  temperature  is  lower  the  lower  the 
input  power  (smaller  p  and  I)  the  lower  the  temperature  of  the 
cooling  jacket,  the  lower  the  molar  fraction  of  Ng  and  the  molar 
fraction  of  the  higher. 


11 


•*  T//77 


1 


Figure  6.  Difference  between  temper-^ 
ature  at  the  center  of  the  tube  and 
at  the  wall,  as  a  function  of  total 
pressure  in  the  tube,  for  various  dis¬ 
charge  currents,  calculated  by  the 
two  methods . 
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Figure  8.  Current  dependence 
of  the  temperature  of  the 
gas  mixture  at  the  center  of 
the  tube. 
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Figure  10  Temperature  of  the 


Figure  7.  Temperature  of  the  gas 
mixture  at  the  center  of  the  tube 
as  a  function  of  total  pressure, 
composition  of  the  mixture  as  para- 
mete'"  ( f r ac t ion  o f  He  )  . 


Figure  9.  Temperature  of  gas  mix¬ 
ture  at  the  center  of  the  tube  as 
a  function  of  the  capillary  wail 


gas  mixture  at  the  center  of  the 
tube  as  a  function  of  Np  content 
for  the  mixture  <;o.  v.  n*  i 
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“  Figure  11.  Temperature  of  the  gas 
-)  mixture  at  the  center  of  the  tube 
J  as  a  function  of  the  content  for  the 
:  mixture  ox :  n,  ;  ne  -- x  .•  i, 


2.2.  Electron  Density 


Experimental  studies  [10]  have  proven  that  the  electron 
density  is  almost  linear  in  the  discharge  current,  the  propor¬ 
tionality  constant  being  a  function  of  the  gas  mixture.  The 
proportionality  constant  is  estimated  using  the  results  in 
[12],  which  show  that  the  drift  velocity  of  the  electrons, 

,  for  Jiliy  in  the  range  (2-5)  x  10”^  V/cm2,  is  in  the 
range  (4-3)  x  10  cm/s  and  depends  only  slightly  on  the  gas 
mixture The  E.N  ratio  is  determined  by  the  necessity  to  keep 
the  ionization  at  a  high  enough  rate  to  balance  the  loss  of 
charged  particles  in  a  self-maintained  discharge,  and  to  opti¬ 
mize  the  efficiency  of  the  vibrational  excitation  in  the  laser 
discharge.  Experimental  measurements  on  a  laser  with  wave¬ 
guide  [17]  give  a  ratio  E/N  of  2.3  x  10 lo  Vcm2  (from  the  slope 
of  the  voltage-pressure  curve).  For  this  value  calculations 
in  [12]  give  r„  ~  i  •  loVm/s  .  If  we  substitute  the  drift 
velocity  in 


where  N  is  the  electron  density,  and  e  the  electron  charge, 

we  obtain  •-=  l,'»- lo*  .  The  optimum  population  inver 

r^j 

sion^is  in  a  narrow  range  of  the  fractional  ionization  Ne/N  = 
10'a  -  lO'^.^As  has  been  previously  shown,  the  electron  distri 
bution  across  the  discharge  tube  is  not  uniform  and  thus  the 
current  density  can  not  be  obtained  from  the  current/area  ratio 
Assuming  the  same  radial  electron  distribution  as  in  the  pre¬ 
vious  paragraph  (Bessel)  and  using  the  integration  of  this  dis¬ 
tribution  from  (6),  one  obtains: 


j  (in. -I /fin4)  -  7h,7h  /(mA)At2(iiini4) 
and  thus  .V,  (cnr:i)  l,lO(j-  lou  /(m.l)/«4(ium2). 

2.3.  Rates  of  Vibrational  Excitation 


Pumping  rates  electron  molecule  for  each  energy  level,  R, 


are  a  function  of  the  transverse  cross  section  of  vibrational 
excitation  by  the  electron,  a  ,  and  of  the  electron  distribution 
function,  which  in  turn,  is  a  function  of  all  the  important 
collisional  processes. 

A*  (17) 

Thus  the  excitation  rates  for  each  vibrational  level  are 
different  for  each  gas  mixture  and  for  each  value  of  E/N.  More¬ 
over  in  the  case  of  N2  molecules,  the  first  3  vibrational  levels 
are  excited.  The  problem  is  simplified  significantly,  if  an 
effective  rate  for  the  vibrational  excitation  by  electrons 
of  the  first  vibrational  level  is  defined. 

Nighan  [19]  calculated  these  rates  for  N2  and  CC>2  as  a  func¬ 
tion  of  the  average  electron  energy  (which  for  a  Maxwell  dis¬ 
tribution  represents  the  electron  temperature).  For  average 
energies  larger  than  1.5  eV,  typical  for  CC>2  lasers,  one  ob¬ 
tains  R.  *  4  x  10"9  Ne,  R„  =  6  xlO-9  N  and  R,  =  2  x  10“^  N, 

1  i  2  o  e  3  s 

where  R  is  in  sec-  and  N  is  in  cm-3. 

e 

2.4.  Rates  of  Collisional  Relaxation 

The  collisional  relaxation  rate  k(s_1)  in  a  CC>2,  N2  and 
He  mixture  is  obtained  from 

4  —  kLO-pt  1 1,  4-  kS*pNt  +  kH'p He) 

or 

*  -  (k^ co,  +  k^St  +  k»'^H.)p, 

where  k1  are  the  collisional  relaxation  rates  in  the  pure 
gases  expressed  in  s-^  torr-1. 

A  survey  of  the  experimental  measurements  of  the  rates  k* 
and  their  temperature  dependence,  is  given  in  [8],  [9].  In  cal¬ 
culating  population  and  gain  the  following  values  for  the  col- 


(16) 

(19) 
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lisional  relaxation  rates  are  used: 


y«uVl)  =  4,197  •  10 


xp  [10  32’(8,84  •  10~*2’  —  2,07)] 


valid  in  the  range  30Q-3000K;  at  300K,  k^  ~  13  ®23  s-1  ton-*1 ; 
All'--  =  t*xp(ll,83  -  7,495  7’-  ‘/a  _  032,14  2’-"*  +  2  23 9T~l), 

co? 

valid  in' the  range  250-2500K;  at  30QK,  =  329  s"1  toir'1 ; 

A*.  -  exp(  — 11,184  -t-  507,831-"*  -  4  574,91-"*  +  12  724T*1)* 

valid  in  the  range  3-00-1000K;  at  300  K,  ktjf  =  ;,‘1  torr-1 ; 

/fsi"  =  txp(24, 538  -  205,482’-"*  +  763,322’-"*  +  718, 832’~l), 


valid  in  the  range  300  -1000K,  at  300K,  k.,,  =  73  s.-'torr-1 ; 

Ay,  =  fcaoxp  (-  hA'/AT)  =  fc.,,exp(- 1382,8  T-1) ; 

A];1*  -=  3,09  10*  2’-*  exp  (-  722’~"*), 


valid  in  the  range  300  -  3000K,  at  300K,  kl02  *  219  s"ltorr_* » 

^ i’ll1  -  4,538-  I0s7’-1(»xp  (— 722*  ,/3) 

N? 

valid  in  the  range  300  -  3000K,  at  300  K,  k^Q  =  32  s_l  torr  1 ;  (27) 
/«•!![  =  4,035-  L0-*  T~x  exp  ( -40,03  J* "*), 

He 

valid  in  the  range  30G-30COK,  at  300K,  k^Q  ~  3  •'*“•**  «’*  torr"1 ;  (28) 

'■-■i  =  A,,  exp  (_<joo,  38 21'1). 

2.5.  Rates  of  Electron  Relaxation 


The  only  levels  which  relax  through  electron  collisions  are 
00°1  in  CC>2  and  v  *  1  in  N^. 

The  relaxation  of  the  C 0^  molecules  from  the  00°1  level 
through  electron  collisions  takes  place  on  one  of  the  levels 
in  group  1  and  not  directly  to  ground  level,  as  was  assumed  in 
[2]  -  [51.  This  rate  is  [20  ] :  =  1,77  •  LO-*  cm3  sr*.  >  thus  for 


iV,  =  3  •  ioncm-3  typical  for  a  waveguide  laser  —  f>,;$  •  iO*s'1 

This  rate  compares  well  with  the  measurements  by  Gower  and 
Carswell  [21]  (r»  •  10s  s"1)  ,  and  is  somewhat  smaller  than  that 
of  Moore  et  al  [1]  cMO*sr‘) 

The  v  =  1  level  of  N2  is  assumed  to  relax  through  elec¬ 
tron  collisions  at  a  rate  =  2  •  io_*  uni3  tr1  •  This 

assumption  was  used  by  Moore  [1]  and  Gordietz  [2]-[4], 

2.6.  Relaxation  rates  through  diffusion  and  collisions 
at  the  walls. 

The  relaxation  rates  through  diffusion  and  collisions  at 
the  wall  in  calculated  from  [3] 

*••„  =  !*  *£>/«*  (29) 

where  D  is  the  diffusion  coefficient,  and'!*'  is  the  solution 
of 


where  o  is  the  average  velocity,  }0(i*)  and  ^(n)  are  the 
Bessel  functions  of  the  first  species  of  order  0  and  1  respect¬ 
ively,  and  r  is  the  reflection  coefficient  of  the  wall  (re¬ 
flection  coefficient  =  1  -  accomodation  coefficient). 

For  a  hollow  circular  wave  guide  equation  (30)  can  be 
simplified  by  a  Taylor  series  development  around  !*  *  2.405 

[7] 

it  =  2,405  SL=J&  ~  2(1  +  r)D 

(l-»r)ra  '  (31) 

The  reflection  coefficient  of  the  wall  was  measured  by 
Kovacs  et  al  [22]  and  is  r  *  0.78.  The  average  velocity  of  the 
CC>2  molecules  depends  on  temperature  only  tf(cm/H)=2,iy5 -iO*}fT(K) 


In  a  C02,  N  , 


He  mixture  the  diffusion  coefficient  is 


obtained  from  [8] 


A  =  V  AV‘  _ 

D  f  «<  ^ 


■  i“)  ‘p“ 


A»2_? oo  \«*  ^ 

A  l  T  J  Aa,’ 


where  is  the  density  of  species  i ,  is  the  collision 
cross  section  and  *  is  a  constant.  For  the  case  of  C02 
and  ,  .  .  \  /  *uu\  \  a/2 


1 

-  f  W*>.  .  _ 

dv, 

VHe  V 

J  3  0  0_Y 

Deo,  " 

=  [  -r* 

•  Deo, -eo. 

Deo,  -  N, 

Deo, -He  / 

A  r  J 

1 

f  'Wo.  _ 

'V,  . 

Z'  300  V 

t3 

1 

V  Ds.-eo, 

Dn.-n, 

DN,-He  i 

{  T  ) 

2 

(3*0 

(35) 


The  constants  are  those  of  dordietz  et  al  [3]  at  />  —  itorr: 

/X’o.-coi  =  !il  tim‘/11  (Fordietz  gives  sect u?/h  ,  but  this  was  corrected 

to  agree  with  the  value  of  0.066  cm  /s  reported  by  Kovacs 
et  al  [22]  at  1  atm  for  CO,  molecules  in  the  00°1  state), 

Deo,  n,  —  ll'O  cm2/a.  Deo,  -  ne  —  .100  cm2/s,  Z>s, -to,  =  ISO  Dn,-n,  =  DM  cm-/s 

/In,  -  He  =  ')t>0  cm*/*' 

Thus  ,;>f  =  ’j»  =  *a>,  W  ( 36 ) 


2.7  Fractions  of  excited  molecules  on  the  upper  and  lower  laser 
levels 


The  fractions  of  molecules  fz  and  on  the  upper  and 
lower  laser  level  are  calculated  from  [233 


m  -  >.  ( -  £;)[£*', «p(-  ^)] ", 


(36) 
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t 


where  f(i)  is  the  fraction  of  molecules  on  the  i-th  vibrational 
level,  expressed  for  a  simple  harmonic  oscillator;  g  are  the 
degeneracies ,  which  for  the  symmetric  and  asymmetric  elongation 
mode  of  C 0,  are  unity  and  for  the  deformation  mode  of  C02, 

(  iyi>  fa  -  q  ~  i  J  is  the  vibrational  temperature. 


In  particular  for  the  00°1  level 


h  = 


>-\p(  -  \Et:kTt.) _ =  l  -evnf  -  - 

lip,'- Isjt r, s7|  I  -  M )  -  Pl  >rJ 


=  l  -  r,, 


(39) 


where 


/  AA\.  \ 
'*  ‘‘XP(  kT,J 


a  a'..  \  .v2  yvib 


-V, 


(40) 


Qvb  is  the  vibrational  partition  function. 

Because  A’./.V,  <  i  .is  a  good  approximation.  The 

/,  a;  1  modes  are  characterized  by  the^ame  vibrational  tempera¬ 
ture,  Tvl.  The  lower  laser  level  relaxes  rapidly  collisionally 
to  the  01^*  level,  half  the  level  (lo°o,  02°0),  .  Thus  the  frac¬ 

tion  f,  on  the  lower  level  is  twice  the  fraction  on  the  level 
0110 


/.  - 


3 


I  (i  +  1)4 


(41) 


where 


(42) 


Summing  in  the  denominator 

f  =  3j'»(l  -  -fi)2 
■'*  3  —  2x, 

for  ru  -  300  K,  x,  =  0,04073  f,  -  0,0387 

f.  has  a  maximum  for  x,  —  0,40092  (Tc,  =  1067,  87 lv), 
/,  =  0,19(i41  .  At  r„  =.  400  K,/,  0,09. 


(43) 

The  function 
,  for  which 


k 
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2.8.  Induced  Emission  and  Absorption 


The  rate  of  induced  emission  in  the  presence  of  monochro¬ 
matic  radiation  with  intensity  /v(W/m*)  at  the  center  of 

the  line,  is  obtained  from  Einstein's  treatment  [24] 


SttA 


'• J 


(44) 


where  A21  is  the  Einstein  coefficient  of  spontaneous  emission 
(.L,  -=  o,j  S-1  for  the  P(20)  line  of  transition  o<ri  —  (ioJo,  02°0M 

.  at(0)  is  x;he  form  function  at  the  center  of  the  line,  and 
}'»  ~  •rjiy-i  is  the  fraction  of  molecules  on  the  upper  laser 

level  (00°1)  which  are  in  the  state  characterized  by  the  quan¬ 
tum  rotational  number  ./  f-  i  is  obtained  from 


where  E,  =  b\-J)he  is  the  energy  of  the  rotational  level  J” 
with  respect  to  the  ground  rotational  level,  and  *V)  - 
BJ(J  -i-l);  Qt*  —  kTrlnteU  is  the  rotational  partition  function,  B 
is  the  rotation  constant  of  the  corresponding  vibrations 
level  and  ?  is  the  molecular  symmetry  number  (mocjuj_e).  2 

Tr  the  rotational  temperature,  assumed  equal  to  the  kinetic 
temperature  of  the  gas.  For  J*  *  20  and  T  s  400  K,  F2  = 
0.063989. 


For  the  C02  laser,  with  wave  guide,  9(0)  -  -K,  where 
v<-  is  the  collision  frequency,  thus  (44)  becomes 


JKa  =  1,03  •  10  IJ' ,/Vl. 


(46) 


To  calculate  the  collision  frequency  for  the  P(20)  line 
one  uses  [25] 


/  rn 

l  r  J  ■Vco*  "**  "1 


(47) 
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where 


'  =  2.2«  •  10,K‘  ^  W>.  =  0,738  aa  =  0,641. 


The  rate  of  induced  absorption  is: correlated  to  that 
of  induced  emission  by 


iV'2  =  (9*/8i)W 2i  =  9,32-10“  7v*yv,. 


(48) 


The  net  exchange  of  molecules  between  the  upper  and  lower 
laser  levels  due  to  induced  emission  and  absorption  is 


W  -  -  "Wi-V,  =  -  A /  v 


(49) 


3.  Equation  system  and  its  solution 


Returning  to  figure  1  we  can  write  the  equation  system 

which  characterizes  the  population  dynamics  for  the  four 

levels  of  our  model 

-  /(*, „VCI  4-  r(k2lX.  —  kl2Xt)  4-  s(k0lN9  —  — 


cU 

-  +  k2,XeN i  +  w, 

u  *v-  =  R„Xa  -  k,'X,N 3  -  4;ixV2  +  + 

d< 

-  k3iX 0^3  ~  fct^^VS  -  -  W, 

<I.V3 
<1/ 


-  i?3*V§  -  kyX{N3  -  lc32X0X3  +  -  k3rN3. 


V  —  V  —  V  --  V 

V  Vn  V 

,l«  —  *’  i  J- 


(.50) 

(51) 

(.52) 

(53) 

(54) 


Coefficient  r  in  (50)  takes  into  account  that  the  col- 
lisional  relaxation  of  C02  molecules  in  group  2  (with  v 3 , 
excitation  results  in  exciting  more  than  one  quantum  of  the 
C02  molecule  from  group  1  (with  vi  and  v*.  quanta  excited). 
Thus  the  reactions  characterized  by  the  global  rates  k^ 
and  k^2  are  [8] 


(55) 


C(X(0O:l)  -i-  J /  (  0.(1110,  03*0)!  -j-  M  +  < 

(  Oa(00el)  -j-  M  —  00.(11*0,  03*0),!  +  i 


(56) 


<  ^..(Oon)  -r  .1/  CO.(20°0  . . .  <H°0)m  +  21  -r  199  cm'1,  . 

\  0  I  ) 

where  M  can  be  C02,  N2  or  He.  Reactions  (55),  (56)  are  char¬ 
acterized  by  r  =  3/2  and  reaction  (57)  by  a  coefficient  r  » 

2  [  5  J  - 


Coefficient  s  in  (50)  taken  into  account  that  the  colli- 
sional  relaxation  of  molecules  from  group  1  takes  place  mostly 
through  reaction  [8],  [9]. 

CO^Ol'O)  M  —  CONOCO)  4-  M  4-  667  cm*1.  '55) 


Because  the  lower  laser  level  has  an  energy  approximately 

twice  that  of  level  01*0,  reaction  (58)  is  equivalent  to  relax- 

v 

ation  of  only  1/2  quantum  /i;  thus  s=  1/2. 


For  the  system  (50)— (5*0  we  have  to  take  into  account 
the  radial  distribution  of  the  electrons  and  thus  the  radial 
Bessel  distribution  of  Mg ,  R1,  R2 ,  R^ ,  N^,  Np  and  N,.  Inte¬ 
grating  these  equations  across  the  section  (assumed  circular) 
of  the  discharge  tube,  at  steady  state  (d/dt  =  0)  leads  to 
a  new  system  of  equations: 


where 


A,  A,/,  •  rl;.,lS1Cl  2  A  i-  v/.0|.\  „  —  A  ,f,  - 
-  O-V/1,  •  T.  irr,  -  o. 

'/-V,  /.,  A,  x.r.  -  I:.,x,rl  ■  /.i2a \r,  ~ 

•  OO'OV',  --  l.-.£lx,x$rt  -  L,,XA:,  -  W<\  -  0, 
b.,x,xj'2  -  kMx,xa(\ 

T  -  l;,l.X3rl  0, 

-'■»  -  -V„  v,r,  r  X.(\, 

v*  -  -V;;  x,rlf 

r:,,-  (-.19V '«),[,•  i  0.4,'H, 

ffn>  l  'l?5uO;C-M0.-,  r/«i)dr  a  0,270, 


r, 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 
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and  N^,  N2  and  in  the  new  system  of  equations  are  the 
axial  densities  (peak). 


Introducing  the  notation  ,ks>j  —  Ai,  =  A si Ci,  k* 
ii,  —  It t(  1,  J-2  =  R..l  u  R3  —  RJ\,  k-ii  —  Ay.*1,  and 

A,  =  A)  -u  A'i.-Ye  -r  Aa„. 

A'.  -  A*.)  -f-  /4.V,  -I-  Ay,  Ay.,, 

_  A1  y-  *'Ay,  -r-  .v/.'ui  -t-  aAlu  -r  Ay,  -r  VV^s/i 
L;1  “  rk  ~~-  lt\  -  #*„,  -  Ay.Y,  -i-  »ra/4 

A' 1  ^A‘m _ 

1  »a-.,,  —  a;  —  ■,<AU1  —  k'-,yt  --  iv  2, /a 

Vtj  -•  AJj(  li'i  —  A'12  -j-  A'jH  j.,  IV  12/i  t  AI0A3), 

A',  •-=  A-3i(A'is  -  -  tC3WiJs  +  11  Vi  -  AV«-'j)  t 

+  /4(A :,iv'21/j  -  a2a'4  -  a,), 

A',  •—  —  A\,2U3  —  k33KlK3, 

Ag  —  A1(A'ij  —  i?2  —  A'3u  jjJj  -  VI  jj/j  K.,K3), 

A,  A-as(Ws"-  A'4U'lt/3  A, A,). 


A'jfA  4  1  » 


(66) 

(67) 

(66) 

(69 ) 


470) 


(71) 

(72) 

(73) 

(74) 


A  A,(  A’j  A.VV  ,,/,  A.  A';). 


(75) 


A,,  /.:l,(  A,  A,). 


(76) 


we  can  determine  the  population  of  groups  1,  2  and  3. 


■v, 


(A  ^  A/  A  -  A  31  *-  A,) 
•JA, 


~  1  (An A',  j-  A 7 A n,  -4-  A'J~  -  4AS( A,,A’,~  1  A'^A'i  A'tl' 

2  A,  '  ‘ 


(77) 


(we  choose  the  solution  with  minus  before  the  square  root, 
because  the  one  with  plus  gives  a  density  of  molecules  in 
group  1  larger  than  the  total  density  of  C02  molecules) 


Ar.  A3A\  +  A',  A',, 

A' 3  -  (Ay.  -  /(■)  -r  W.JiVA'jA',  -  (A.  ‘  WtJi)Xs!Kx  - 

r  (/fj.Vji  ►  AjAij/v,. 


(78) 

(79) 
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and  the  number  of  CO ^  and  molecules  in  the  ground  state 
using  equations  (62)  and  (63) 

'Thus  the  solution  of  the  system  of  equations  which  charac¬ 
terize  a  CO2  laser  with  wave  guide  allow  the  determination  of 
the  population  inversion  and  the  population  of  the  laser  levels. 


The  small  signal  gain  at  the  center  of  a  line  in  one  P 
branch  is  calculated  from  these  populations  (witn  /, 
and  H'l#  =  W.,  —  in  [24] 

K'F(v.J)J  'J  P?i~sex|»f -J(J  -  l)  -  t.xp|  ■•/(./  ■  1)  ML 

''I)  L  t.  ‘r\\ 


where 


i  l»TJ  V„ 

!/,  tv,; 


-N  1  > 


(31) 


where  K\  =  ■  U)~-D  is  the  matrix  element  of  the  circle 

moment  for  vibrational  transitions ,  //,  =  o.duoik  cm**,  h,  ==  <>.;S8714  cm*1, 
vo  =  2,480881  •  iolall2  is  the  central  frequency  of  the  bar..!  ,  *’,r’  ^ 

—  1  -<-o,OtH)ik/  —  OJHJVOH'/3  is  the  coefficient  of  interaction 
vibration  rotation,  v,  -=  •  lo,JIIz  is  the  frequency 

of  line  P(20),  and  iar  0  --==  he  IS  Ik  is  the  rotational  tempera¬ 

ture  [26]  -  [23],  Introducing  these  constants  in  (30)  and  (3l) 
for  line  P(20)  of  transition  00°1  -  (10°0,  02°0)I  we  obtain 


/<l/n  i,&2 2-10-*  ,  0-5 

;  = — — T !-N i/>M>(-2ll,7/rv  l,0078.V,/lo.x1»(-l*;r,,8/7'))  (02. 

\  ■"  !  v,  l 


In  the  exact  calculation  of  the  small  signal  gain  one 
has  to  take  Into  account  the  losses  in  the  wave  guide..  The 
losses  in  the  propagation  mode  EH^  in  the  cylindric  dual 
wave  guide  is  given  by  [29] 


*11 


Kf(v„), 

«a 


(63) 


23 


where  "m  —  is  the  first  root  of  equation  •]„(«.*)  =u, 

>=io,59  ,um  is  the  wavelength  of  the  radiation  propagating 

through  the  waveguide  and  v„  (■/■  u-  DiC2]f '/•  ^T).  ,  where  / 

is  the  complex  index  of  the  wall  material  ( v  =  l  .910 -4- 0.077 
for  Pyrex  glass  [30]),  thus  [£e(v„)-=l,id«  .  Under  these  con¬ 

ditions  5 tn(dli/m)  =  0,105/«3(mm3) 

4.  Results 

Populations  of  the  laser  levels,  small  signal  gain  and 
saturation  intensity  were  calculated  for  various  regimes  of  the 
CO^  laser  with  wave  guide.  Most  of  the  calculations  were  done 
for  the  mixture  ro. :  N. :  1U>  =  1 : 1 :  ;j,  characteristic  for 

lasers  with  gas  circulation  [31],  but  other  mixtures  were  tested. 

Figure  12  shows  the  variation  of  the  population  -V]  . 

v. 

as  a  function  of  the  discharge  current  for  three  different 
pressures,  for  a  1  mm  diameter  waveguide.  For  100  and  150 
torr  pressure,  the  population  inversion  is  maintained  over 
the  whole  range,  for  30  torr  the  population  inversion  occurs 
only  for  currents  smaller  than  5  mA. 

The  pressure  dependence  of  N2/N^  for  various  discharge 
currents  is  shown  in  figure  13. 

For  a  1  mm  capillary  and  CO. :  :  Uo  =  1 : 1  s  3,  the 

maximum  of  the  N2/N^  ratio  occurs  in  the  pressure  range  100 
to  150  torr. 

To  verify  our  model  the  small  signal  gain  was  calculated 
for  the  experimental  conditions  used  by  Abrams  and  Bridges  [32], 
this  comparison  is  presented  in  figure  14.  The  agreement  is 
very  good  as  far  as  the  absolute  value  of  the  gain  is  concerned, 
the  experimental  curves  are  shifted  toward  higher  pressures 
though.  This  shift  can  be  explained  by  the  use  of  berillia 
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Figure  12.  Absolute  populations 
of  groups  1  and  2  as  a  function 
of  the  discharge  current  in  a  1 
mm  capillary. 

\ 


capillaries  in  [32]  for  which  the  cooling  is  more  efficient 
and  by  relaxation  rates  at  the  wall  different  than  those  used 
in  calculation  (for  Pyrex  glass). 

Figures  15  and  16  present  the  small  signal  gain  as  a  func¬ 
tion  of  pressure  (gas  mixture  as  parameter)  and  the  radius  of 
the  wave  guide  respectively.  In  figure  15  the  gain  was  calculat¬ 
ed  for  mixtures  in  which  the  fraction  of  He  is  modified. 

If  for  the  1:1:0  mixture  the  gain  decreases  monotonously  and 
rapidly  with  pressure,  the  mixtures  containing  He  show  a  maxi- 
mu.  The  agreement  between  experimental  measurements  of  the 
gain  for  a  laser  with  a  =  1  mm,  I  ^  s  mA  and  : 

Hi*  =  L :  1  :  [17]  and  our  calculations  is  very  good. From 

figure  lo  in  which  the  dashed  lines  represent  the  calculated 
gain  when  neglecting  losses  in  the  waveguide,  and  the  continu¬ 
ous  one  which  takes  losses  into  account,  we  concluded  that 
the  optimum  diameter  for  the  capillaries  is  in  the  range  1-1.5 
mm  and  that  losses  are  important  only  for  diameters  smaller 
than  1.5  mm. 

Aside  from  calculating  populations  and  gain  the  proposed 
model  can  be  used  to  calculate  the  saturation  intensity  (satur¬ 
ation  parameter).  Assuming  that  in  waveguide  lasers  the  line 
widens  homogeneously  with  pressure,  the  saturation  intensity 
is  determined  as  the  value  of  intensity  for  which  the  small 
signal  gain  decreases  to  half. 

Figure  17  shows  the  pressure  dependence  of  the  saturation 
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Figure  13.  Population 
ratio  for  groups  1  and  2 
as  a  function  of  total  pres¬ 
sure  (a  =  0.5  mm). 
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Figure  14.  Comparison  between 
calculations  of  gain  and  experi¬ 
mental  measurements  of  Abrams 
and  Bridges  [32]  for  the  same 
discharge  conditions. 
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Figure  15.  Small  signal 
gain  as  a  function  of 
total  pressure  in  a  tube 
of  2  mm  diameter. 
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dependence  o/*the  small  signax  gain 
for  various  mixtures  and  discharge 
currents . 


intensity  for  various  discharge  currents.  Saturation  intensity 
increases  with  pressure-  which  is  also  shown  by  simpler  models. 
The  large  values  of  parameter  I_  for  waveguide  lasers  can  be 

o 

explained  by  the  large  pressures  used  in  these  lasers,  and  be¬ 
cause  the  process  of  diffusion  and  relaxation  at  the  wall  is 
more  efficient  than  for  conventional  lasers  (low  pressures). 


The  dependence  of  the  current  on  saturation  intensity  for  a 
waveguide  2  ram  in  diameter  is  shown  in  Figure  18.  Agreement 
with  the  simpler  models,  which  lead  to  linear  increases  in  the 
parameters  of  saturation  with  current,  is  likewise  good. 
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Figure  17.  Pressure  dependence  of  the 
saturation  intensity,  discharge  current 
as  a  parameter. 


Figure  18.  Discharge  current 
dependence  of  the  saturation 
intensity  for  a  2  mm  tube. 


5.  Conclusions 


A  four  level  system  is  used  to  model  the  physical  processes 
and  equations  characterizing  the  C02  lasers  with  waveguide.  The 
system  of  equation  is  solved  by  computer.  It  allows  to  deter¬ 
mine  the  dependence  of  important  parameters  such  as:  kinetic 
temperature,  population  of  laser  levels,  small  signal  gain, 
saturation  intensity, as  a  function  of:  discharge  current,  total 
pressure,  gas  mixture  etc.  The  agreement  between  calculations 
and  experiment  is  satisfactory. 

These  calculations  can  be  used  to  optimize  the  functioning 
of  C02  lasers  with  waveguide  according  to  their  use  (small  size, 
efficient  cooling,  high  oscillation  band)  .'  ^f Though  the  model 
is  applied  to  waveguide  lasers  it  can  be  extended  to  other  types 
of  CO,,  lasers. 
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